Neutron stars in Einstein-aether theory 



Christopher ElingD, Ted Jacobsonlj, and M. Coleman Miller@ 



* Department of Physics 
University of Maryland 
College Park, MD 20742-4111, USA 

^Department of Astronomy 
University of Maryland 
College Park, MD 20742-4111, USA 



Abstract 

As current and future experiments probe strong gravitational regimes 
around neutron stars and black holes, it is desirable to have theoret- 
ically sound alternatives to General Relativity against which to test 
observations. Here we study the consequences of one such general- 
ization, Einstein-aether theory, for the properties of non-rotating neu- 
tron stars. This theory has a parameter range that satisfies all current 
weak-field tests. We find that within this range it leads to lower max- 
imum neutron star masses, as well as larger surface redshifts at a par- 
ticular mass, for a given nuclear equation of state. For non-rotating 
black holes and neutron stars, the innermost stable circular orbit is 
only slightly modified in this theory. 
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1 Introduction 



General Relativity (GR) has passed every test so far, from solar system dy- 
namics and light bending to precision measurements of the orbits of binary 
pulsars PQ. Nonetheless, tests in strong gravity remain elusive, largely be- 
cause phenomena in strong gravity involve other uncertain physics as well. 
For example, spectral profiles of Fe Ka fluorescence lines in active galac- 
tic nuclei and stellar-mass black holes [21 111 El E] are consistent with the 
expectations of gas streamlines near rapidly rotating black holes. However, 
precision tests are not yet possible because of unknowns about the emission 
profile and other complications. 

More robust tests are on the horizon. Larger-area X-ray detectors such 
as Constellation- X [7] may be able to track the motion of individual emitting 
elements in a disk, hence mapping out the spacetime near a rotating black 
hole. Even more robustness is likely to come from detections of gravitational 
waves from double black hole mergers of various masses and mass ratios seen 
with ground-based interferometers such as LIGO [5], VIRGO [5], GEO-600 
[lOj . and TAMA [11] and later with space-based instruments such as LISA 
|12j . A third source of strong field tests could exploit neutron star masses 
and surface redshifts. It is that sort of system that is the subject of this 
paper. 

It is therefore desirable for current and future data to have theoretically 
sound — and preferably well-motivated — alternatives to or generalizations of 
GR. The simplest alternatives that have been considered are various scalar- 
tensor theories. In the case of the well-known Jordan-Brans-Dicke theory 
it has been shown [13J that the predictions of the theory in both the weak 
and strong field regimes deviate from GR by a parameter that is tightly con- 
strained by post-Newtonian Solar System experiments. Thus, the properties 
of compact objects such as neutron stars in Jordan-Brans-Dicke theories that 
pass these weak field tests must be very close to those found in GR. However, 
Damour and Esposito-Farese [LT] found a wide class of other scalar-tensor 
theories exhibiting "spontaneous scalarization," where weak field constraints 
are met, but compact objects have significant deviations from GR in the 
strong field regime. Recently, [15] studied the properties of non-rotating 
neutron stars in these theories, finding larger stellar masses than in GR and 
larger surface redshifts for a given equation of state. These results were 
then used to put an observational constraint on one of the parameters of the 
model. 
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Over the past several years, a number of further alternatives to GR have 
been proposed in the context of incorporating Lorentz violation (LV) into 
gravity; see for example Refs. [HI HZl UHl UHl I2Q1 [2H [22l [23] and references 
therein. These theories can be thought of as low energy effective field theory 
descriptions of local Lorentz violation (LV) possibly arising from quantum 
gravitational physics at energies near the Planck scale. In this paper we study 
one of these models, "Einstein-aether" theory (or "ae-theory" for short) [24], 
in which a dynamical unit timelike vector field u a is coupled to gravity. The 
LV here is a sort of spontaneous symmetry breaking. The action is Lorentz 
invariant, but in any solution defines the 4- velocity of a local "preferred" 
frame at each spacetime point, which breaks local boost symmetry. 

The general ae-theory dynamics is dependent on four dimensionless cou- 
pling parameters. Weak field constraints on ae-theory can all be met by 
restricting these parameters to a subset of the four dimensional parameter 
space that is very narrow compared to unity in all but one dimension. The 
nature of these constraints is summarized in the following section. They are 
derived in the weak field regime, except for strongly self-gravitating neutron 
star binaries. It is therefore interesting to ask what are the deviations from 
general relativity in the strong field regime, and to compare that behavior 
with observations of astrophysical systems. 

The main objective of this paper is to study the properties of non-rotating 
neutron stars in ae-theory. The stellar solutions depend on only one com- 
bination of the theory's coupling parameters. We consider six candidate 
equations of state, three with purely nucleonic degrees of freedom and differ- 
ent hardness, and three involving quark matter with different bag constants. 
By numerical solution of the field equations interior to the star we obtain the 
maximum mass, relation between mass and radius, and surface redshifts, all 
function of the coupling parameter. 

It turns out that the maximum neutron star mass is less than in the case 
of GR, and is smaller for larger values of the coupling parameter. Thus, once 
the equation of state becomes known well enough, it will be possible to place 
an upper bound on the coupling parameter by observations of neutron star 
masses. Nonstandard relations between mass and surface redshift also occur, 
providing another possibility for interesting phenomenology and constraining 
the coupling parameter. 

We also examine the location of the innermost stable circular orbit (ISCO) 
as a function of mass, to determine whether that might provide further useful 
observables distinguishing GR and ae-theory. We find however that the ISCO 
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is nearly unchanged for reasonable coupling parameters and non-spinning 
objects. 

The structure of this paper is as follows. In Section II we first review the 
basics of Einstein-aether theory and summarize the current constraints on 
theory. In Section III the equations of structure for fluid stars in ae-theory 
are presented, together with the form of the analytic exterior solution to 
which the interior must be matched. Using the exterior solution, expressions 
are obtained for the surface redshift and ISCO radius that can be employed 
with the numerical solutions to obtain the observable quantities. We present 
the numerical results in Section IV, together with the constraints on the 
coupling parameter that can be obtained with these results. In Section V 
we conclude with a brief discussion of prospects for further constraints from 
more precise neutron star measurements and rotating black hole solutions. 

2 Einstein-aether theory 

The action for Einstein-aether theory is the most general generally covariant 
functional of the spacetime metric g ab and aether field u a involving no more 
than two derivatives (not including total derivatives), 



16ttG 

where 



j V~9 (L x + Lmattcr) <&X (1) 

-R - K ahmn W a u m W h Un - \(g ab u a u b - 1) (2) 



ab 
mn 



and ^matter denoted the matter lagrangian. Here R is the Ricci scalar, K, 
is defined as 

Ktn = cig ab g mn + c 2 6 a m 5 b n + c 3 5 a J b m + c A u a u b g mn (3) 

where the Cj are dimensionless coupling constants, and A is a Lagrange mul- 
tiplier enforcing the unit timelike constraint on the aether. The convention 

used here for metric signature is (H ) and the units are chosen so that 

the speed of light defined by the metric g ab is unity. Note that since the 
covariant derivative operator V a involves derivatives of the metric through 
the connection components, and since the unit vector is nowhere vanishing, 
the terms quadratic in Vw also modify the kinetic terms for the metric. 
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The matter Lagrangian L matter generically will be a functional of a collec- 
tion of matter fields (denoted as if)) along with g ab and u a . However, following 
the observational constraints on Lorentz violation in the matter sector, we 
assume here when studying the neutron star solutions that there is no sig- 
nificant coupling of matter to u a . The absence of coupling of u a to matter 
has no theoretical justification in this purely phenomenological approach, 
and may be regarded as unnatural. However our goal here is just to explore 
consequences of gravitational Lorentz violation in a phenomenologically vi- 
able setting. It remains an open question whether this can emerge as an 
approximation to a more fundamental underlying theory. 

The field equations from varying ([T]) with respect to g ab , u a and A are 
given by 

G ab = T^+MGT M ab (4) 

V„ J a m - C4UaV m U a = A« m , (5) 

g ab u a u b = 1, (6) 

where 

J a m = K ab mn V b u n . (7) 
The aether stress tensor is given by 

T {u) ab = V m (J ia m u b) -J m [a u b) -J {ab) u m ) 

+ci {{V m u a ){V m u b ) - {V a u m ){V b u m )} 
+c 4 u a ii b 

+ [u n iy m J mn ) - c A u 2 ] u a u b 
1 

— -^L u g ab , (8) 

where L u = -K ab mn V a u m V b u n . The La grange multiplier A has been elimi- 
nated from ([8]) by solving for it via the contraction of the aether field equation 
© with u a . 

2.1 Observational constraints on the parameters a 

In the weak-field, slow-motion limit ae-theory reduces to Newtonian grav- 
ity (25], with a value of Newton's constant Gn related to the parameter G in 
the action ([T]) by 

G 

GN = l-(c 1 + c 4 )/2- (9) 
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The phenomenology of Einstein-aether theory has been extensively stud- 
ied over the last few years. Theoretical and observational constraints on 
the coupling parameters q have been determined from parameterized post- 
Newtonian analysis [261 123 [28] , stability and linearized energy positivity [291 
|30~1 I3T1 |32j, primordial nucleosynthesis [25], and vacuum Cerenkov radia- 
tion [31J. The combined constraints from all of these are reviewed in Ref. [28J. 

To summarize here, all parameterized post-Newtonian (PPN) parameters 
except the preferred frame parameters ai,2 agree with those of GR for any 
choice of the q. Observations impose strong constraints on a\ (< 10~ 4 ) and 
«2 4x 10~ 7 ). These parameters can be set to zero in Einstein-aether theory 
by imposing two conditions, on the q, which can be solved to determine 

c 2 = {-2c\ - cic 3 + c|)/3ci 

c 4 = -4/ Cl . (10) 

With this choice the gravitational constant appearing in the cosmological 
Friedman equations agrees with that appearing in the force law between iso- 
lated masses, so there is no further nucleosynthesis constraint. The stability, 
positive energy, and vacuum Cerenkov constraints then impose the inequali- 
ties 

< c+ < 1 

< c_ < c+/3(l -c+), (11) 

where c± = C\ ± c 3 . 

Further constraints have been obtained using radiation damping in binary 
pulsar systems [33J. An analysis neglecting strong self-gravitating effects 
found that when fflOl) hold, just one condition A(c\, C3) = 1 makes the lowest 
order radiation rate in ae-theory identical to that of GR. This condition is 
satisfied entirely in the region allowed by flTTjfl However, the neutron star 
sources are strongly self-gravitating. It turns out [31] that as long as q < 0.1, 
the strong field corrections are negligible, but for larger coupling values the 
precise radiation damping constraints are not yet worked out. They will lead 
to a modified condition A'(c\, C3) = 1 that will depend on the nature of the 
compact objects in the binary. 

4 In [33] , the A = 1 curve does not fall entirely in the otherwise allowed region, but this 
is due to an error in the analysis there that has since been corrected [34] . 
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Figure 1: Graphical representation of (1111) (solid line) and (jl2jl (dashed lines). 
The allowed region of the parameter space is below the solid curve, above 
c_ = and to the left of c + = 1. The radiation damping constraint will 
restrict to a nearly one-dimensional subset of this region. 
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Neutron star structure constraints along the lines discussed in the present 
paper should eventually be able to restrict c u = c\ + c 4 , which is given by 

C14 = 2c+c_/(c+ + c_) (12) 

when the PPN equivalence conditions ffTUl) hold. The only previous con- 
straint on C14 was the requirement that it be less than 2 in order to maintain 
positivity of Newton's constant 

Fig. [1] shows the region in the (c + ,c_) parameter space allowed by the 
above constraints (other than the radiation damping constraint), along with 
C14 contours. Note that without a constraint on cu, c_ can grow arbitrarily 
large as c + — > 1. Any upper bound on cu will cut off this region however. 
A cu contour intersects the right hand boundary (c + = 1) of the allowed 
region at c_ = cu/(2 — c u ), and intersects the upper boundary at c_ = 
(4/3)c 14 /(2-c 14 ). 

3 Neutron Stars 

Time independent spherically symmetric solutions in ae-theory were exten- 
sively studied recently in a pair of papers on fluid star [35] and black hole [36] 
configurations. For black holes, the aether vector has a radial component, 
but for the case of a star it does not, i.e. it is aligned with the static frame 
defined by the Killing vector, both inside and outside the star. Thus, unlike 
in GR, the exterior solution for a star is not the same as for a black hole. 

In Ref. [35] the static vacuum exterior solution was found analytically, and 
numerical integration was employed to find interior stellar solutions for the 
test case of a constant mass density equation of state. The mass was found 
by matching the interior to the exterior solutions. Unlike constant density 
stars in General Relativity (GR), the mass does not increase monotonically 
with central pressure; rather, there is a maximum mass at a finite central 
pressure beyond which the stars are unstable. This maximum mass is smaller 
than in the GR case. 

It is adequate for our present purposes to restrict attention to non- 
rotating stars, since the effect of rotation on the maximum mass, surface 
redshift, and ISCO is very small for the observationally relevant spins. For 
example, presuming the fractional change in maximum mass scales with the 
square of the spin, Tables 4 and 5 in Ref. [37] indicate that even for a one 
millisecond period the increase of maximum mass is less than 5% in GR. 
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Barring an unexpected much greater sensitivity to a small spin in ae-theory, 
the results we find here for non-rotating stars should be quite reliable except 
for the fastest spinning stars. 

3.1 Stellar Equations of Structure 

The static, spherically symmetric form for the metric and aether can be 
written in Schwarzschild-like coordinates as 

ds 2 = e A{r) dt 2 - B(r)dr 2 - rW, (13) 
u = e- A(r)/2 d t . (14) 

Note that since the aether is aligned with the timelike Killing vector it is 
completely determined by the metric. In particular, the aether is at rest with 
respect to the static frame at infinity, which means that the star is taken to be 
at rest with respect to the aether. When comparing theory and observation, 
it is typically assumed that the background aether frame coincides with that 
of the cosmological fluid. Any particular star will of course have some proper 
motion with respect to this frame, so strictly speaking the physically relevant 
solutions are not of the form (11311 141) . However, assuming a relative velocity 
of order 10 -3 , this discrepancy should not be significant for comparisons of 
much less precision such as concern us here. 

It was shown in [35] that for configurations of the form (I13|ll4p the C2 and 
C3 terms in the action ([T]) and their variations are zero, and thus they do 
not contribute to the field equations. Also, the effect of the c 4 term can be 
absorbed by the replacement c\ — > c\ + C4. Hence the only coupling relevant 
to these solutions is C14 = c\ + C4. The fluid stress tensor appearing in the 
metric field equation (T4J) is 

TS = W) + P(r))v a v b - P(r)g ab (15) 

where v a = e~ A ^ 2 (dt) a is the fluid 4- velocity, p(r) its mass density, and P(r) 
its pressure. 

The metric field equation and the Bianchi identity together imply that the 
sum of the aether and fluid energy-momentum tensors is divergenceless. In 
addition, since the aether does not couple directly to the fluid, its stress tensor 
is independently divergenceless when its field equation and unit constraint are 
satisfied. Therefore the fluid stress tensor is also independently divergenceless 
in any solution. Thus, an appropriate system of equations for the aether plus 
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fluid case is the (i) metric field equation, (ii) aether field equation, (iii) radial 
component of V a Tjf = 0, which is the hydrostatic equilibrium equation for 
the fluid 

P' + \A\p + P) = 0, (16) 

and (iv) an equation of state p = p(P). 

It turns out that the aether field equation §5§ has only a t component, 
which just determines the Lagrange multiplier A. The tt, rr, and 99 compo- 
nents of the metric field equation then imply respectively 

B 1 B' 
= -1 + B + r— - v(8rA' + r 2 A' 2 - 2r 2 A'— + Ar 2 A") - pr 2 B (17) 
B B 

= 1 - B + rA' + vr 2 A' 2 - Pr 2 B (18) 
= ^(2rA'-2r^ + r 2 A' 2 -r 2 A'^ + 2r 2 A")-vr 2 A ,2 -PB, (19) 



where the symbol 



" = ^ (20) 



is introduced to compactify the notation, and we have adopted units with 
8ttG = 1. 

Equation (fT8|) can be used to solve for B, 

B = (l + r 2 P)-\l + rA' + vr 2 A 12 ). (21) 

After substituting this result into (fTTl) and (fT9l) . A" can be eliminated from 
this pair of equations, yielding an equation involving A, A', P, P', and p(P). 
Due to its complexity it does not seem illuminating to display it here. This 
equation combined with ffTBT) can then be numerically integrated to solve for 
P(r) and A(r) starting with initial values at the origin r = 0. (It is possible to 
eliminate A', leaving one Tolman-Oppenheimer-Volkoff (TOV) type equation 
for P{r). Since this TOV equation is quite complicated and doesn't aid the 
numerical integration procedure we will not display it here.) 

To numerically integrate outward we find the power series solution to the 
equations (1161) and (I17lll8|ll9p in the vicinity of r = 0, which is a singular 
point for the equations. In this solution the central value for the pressure 
P(0) = Po is the only free parameter to be specified (A(0) is arbitrary due 
to scaling freedom of the t coordinate, so can just be set to unity). The 
numerical integration can then be started at a small value of r using the 
power series for initial data, and continued to the value r = R, which is the 
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surface of the star where the pressure and mass density drop to zero. At this 
point the values of A and A' > can be matched to the static vacuum aether 
solution discussed in (35] and summarized in the next section. 



3.2 Vacuum solution 

In the exterior (1211) becomes 

B = 1 + rA' + vr 2 A' 2 , (22) 

and the remaining field equations (|17|) and (fT9"j) can be reduced to the second 
order ordinary differential equation (ODE) 

r 2 A" + 2rA' + r 2 A' 2 + v r 3 A' 3 = 0, (23) 

With the substitution 

Y = rA', (24) 

(123]) becomes a first order equation for Y(r) that can solved, after which (j2I 
can solved for A(Y). In terms of the roots of B, 



the result is 



and 



Y ± = (-1 ± Vl^Iu)/(2u), (25) 



B = u(Y-YJ){Y-Y + ), (26) 
"«e*=(^r* (27) 



(28) 



\*-Y/Y+/ 
Y \ /Y-Y\^ 



Y-Y_ \Y-Y 



where r m ; n is an integration constant. Thus the complete solution is known 
up to the inversion of the function on the right hand side of ( 1281) . 

As in GR, the solutions in this family are all asymptotically flat. The 
limit Y — ► corresponds to spatial infinity, where the limiting form of the 
solution is 

B = 1+Y + ... (29) 
N = 1-Y + ... (30) 
Y = r g /r + ---. (31) 
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Here r g is the gravitational radius, which is related to the gravitational mass 
M appearing in the Newtonian potential by 

r g = 2G N M. (32) 

The relation between r g and the minimum radius^! r m j n is given by 

r min /r 9 = (-Y + )-\-l - r + )(WV(W). (33) 

To leading order in 1/r this solution agrees with the Schwarzschild solution 
of GR. The total gravitational mass M of the fluid star can be read off from 
AH together with ([28]), using the definition ([21]), Y(R) = RA'(R). 

3.2.1 Surface redshift 

The light emitted from the surface of the neutron star is redshifted as it 
climbs away to a distant observer. From (|13|) . the surface redshift factor z is 
given by 

z = [N(R)}- 1/2 - 1, (34) 
which can be evaluated directly from the numerical solution using ff271) and 

(El. 

3.2.2 ISCO 

The orbits in the metric (fT3"j) have conserved energy e = Nt and angular 
momentum £ = r 2 {p, where N = e A and the overdot stands for derivative 
with respect to proper time. Since the parameter is proper time, the four- 
velocity has unit norm. This condition can be expressed in the form 

f 2 = V(r) = B^W, (35) 

with 

W = W{r; e, £) = A^" 1 (r)e 2 - r~ 2 £ 2 + 1. (36) 

The ISCO is determined by the conditions V — V — V" — 0, or equivalently, 
W = W = W" = 0, where the prime stands for derivative with respect to 

5 In the pure vacuum solution, there is a minimal 2-sphere "throat" at r = r m i n . In a 
fluid star solution the radius of the star is always greater than r m i n , so there is no throat 
in the spatial geometry [55] , 
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r. Thus the metric function B plays no role. These equations determine r, e 
and I at the ISCO. After some manipulation of the equations we obtain 



F ISCO = Zl±yi±l (37) 



With this result, the radius of the ISCO can be found from ( 1281) and ( T33T) 
given the mass. Expanding in v = C14/8 we find (in units with G-^M = 1) 

nsco ^ 6(1 + [ln(3/2) - l/6)v) ~ 6(1 + 0.030 c 14 ), (38) 

dropping 0(v 2 ) terms. This linear approximation is extremely accurate: the 
relative error grows monotonically from to only about 0.3% over the entire 
allowed range of cu from to 2. 

The angular frequency of an orbit with respect to time at infinity is given 
by uj = <f/t = (£/e)(N/r 2 ). The circular orbit condition yields N£/e = 
(iVr 3 /2) 1/2 = r(YN/2y/ 2 , so u = ^(Y N /2) 1 ' 2 . Expanding again in u, the 
frequency at the ISCO is found (in units with G^M = 1) to be 

u lsco ~ ^=(1 + [-21n(3/2) + l/2]z/) ~ ^=(1 - 0.039c 14 ), (39) 

dropping 0(v 2 ) terms. 

Thus, even for the maximum value C14 = 2, the location of the ISCO is 
only about 6% larger than its value in GR for a star of the same mass, and 
the orbital frequency is about 8% smaller. Since ae-theory agrees so closely 
with GR on these quantities, it is unlikely that in the near future any useful 
constraints can be obtained from their behavior for slowly rotating stars. 
However, as we discuss in Section [5], for the ISCO of a rapidly rotating black 
hole the deviations from GR may well be much more pronounced. 



4 Numerical results 

Here we will compare the properties of neutron stars in GR and ae-theory 
using three hadronic and three quark equations of state (EOS). We label 
these according to whether they are softer (s), medium (m), or harder (h), 
by 

Hs, Hm, Hh <-> A18, A18^UIX, A18UIX (40) 
Qs, Qm, Qh <-> (90, 0), (60, 200), (60, 0). (41) 
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The hadronic models are discussed in [35], and the quark models are MIT bag 
models [39j determined by two parameters (B,m s ), with the bag constant B 
measured in MeV/fm 3 and the strange quark mass m s in MeV. 

The pressure and mass density data tables for these models [ID] were con- 
verted from cgs units to geometrized units, i.e. replacing the energy density 
p and pressure p by G-^p/c 4 and G^p/c 4 respectively, yielding quantities with 
dimension inverse length squared. A curve fitting procedure was then used to 
generate an equation of state function p(P) suitable for the numerical inte- 
gration. The field equations (I17til9j) are written in units with 8ttG = c = 1, 
so to apply them we first multiply the density and pressure in the above 
geometrized units by SnG /G^ = 8tc(1 — c 14 /2). 

In GR five of the six equations of state have associated M versus Pq 
curves containing a maximum mass extremum and regions of stability and 
instability. The exception is the softest quark EOS, Qs, which appears to 
asymptote to its maximum mass value. The GR maximum mass values 
for the Hm and Hs equations of state we find here (2.20M Q and 1.67M Q 
respectively) agree very well with the results obtained in [38] • The maximum 
mass value in ae-theory grows smaller and occurs at smaller values of central 
pressure as cu is increased. For the Qs EOS extrema begin to develop in the 
M versus Pq curve as cu approaches 1. Fig. [2] shows a plot of M vs. R for 
the Hm EOS, each point being determined by a value of Pq, for several values 
of Cx4. As Pq increases the mass values increase sharply to peaks and then 
gradually fall off. The region of the curves for small Pq and larger R up to 
the mass maximum describe stable equilibrium configurations. Beyond the 
maximum the neutron stars are unstable. In ae-theory the minimum radius 
where the equilibrium configuration is stable decreases as cu increases. 

A plot of the maximum mass values for the six equations of state consid- 
ered in this paper is shown in Fig. [3J Horizontal lines mark the certain lower 
bound of 1.44 M® and a benchmark value of 2.0 M & . The dependence of the 
maximum mass on cu is very close to linear for the quark models, with the 
mass changing by roughly 6% as cu increases from to 1. For the hadronic 
models it is roughly linear but steeper, decreasing by roughly 15% over the 
same range of Ci 4 . 

4.1 Maximum mass constraints 

The most straightforward constraint on ae-theory comes from comparing the 
maximum mass values generated with the six equations of state to observa- 
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Figure 2: Total mass versus R for the Hm equation of state for Pq up to 
100 and cu = 0, 0.05, 0.2, 0.5, 1 beginning with the upper (GR) curve. The 
vertical axis is units of solar masses and the horizontal in km. The GR curve 
reaches its maximum mass of 2.2OM at slightly more than 10 km. As c±4 
increases to 1 the maximum mass decreases and the value of the radius at 
the maximum falls to about 9.5 km. 
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Figure 3: Maximum mass vs. C14 for the six equations of state. The hadronic 
models are plotted with solid lines, while the quark models are dashed. The 
thick solid horizontal lines represent the bare minimum constraint of 1.44 M Q 
and a possible constraint value of 2 M Q . 
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tions of neutron star masses in binary pulsars. These masses are not directly 
measured, but inferred from the timing data from a binary pulsar system. 
This data contains information on the Keplerian and post-Keplerian param- 
eters of the system, which depend on the unknown masses tjia and of the 
neutron stars. A determination of the Keplerian and two post-Keplerian pa- 
rameters, such as the secular rate of periastron advance and the magnitude of 
the Shapiro delay, result in two curves in a (wajTOb) mass plane. The value 
of the two masses is the intersection of these curves. Since we are consider- 
ing the subset of Einstein-aether theories that satisfies the post-Newtonian 
constraints, any corrections to how the post-Keplerian parameters depend 
on the pulsar masses would only appear at higher order and is therefore very 
small. Thus, the masses can be inferred as in GR to a good approxima- 
tion. If the maximum mass predicted by Einstein-aether theory is smaller 
than an observed neutron star mass then the theory is ruled out. Currently 
the largest reliable observed mass value is 1.44, M Q from the PSR 1913+16 
Hulse- Taylor binary system. 

There are, however, suggestions from various data that neutron star 
masses can be at least ~ 2M . The neutron star in Vela X-l has an es- 
timated mass of 1.88 ± 0.13 M & [H], and PSR J0751, a pulsar in a detached 
low-mass binary, has a reported mass of 2.1±0.2 M & [12]. Furthermore, there 
are indications (although not as definitive) for neutron star masses greater 
than or of order two solar masses in several low-mass X-ray binaries based 
on the inference of the orbital frequency at the ISCO from their kilohertz 
quasi-periodic brightness oscillations (QPO's)[l3l HU 113 EH]; for an alter- 
native view, see [UJ. At the < 700 Hz spin frequencies of these stars, the 
dimensionless angular momentum is only 0.1 to 0.3 [37], so this would imply 
that the ISCO is obtained from a near-Schwarzschild spacetime. Since the 
ISCO we find in Section 13.2.21 is very close to the GR value, the derived 
mass should be the same in ae-theory as in GR to a good approximation. 
As a benchmark, we will consider the limits on cu that would result from a 
measured gravitational mass of 2M . 

Fig. [3] shows that in GR (ci 4 = 0) all six equations of state respect the 
lower bound of 1.44 M Q solar masses. For four of the equation of state models 
the 1.44 M & mass cutoff does not yield any constraint on C14. For the Hs and 
Qs EOS there are weak constraints that C14 be less than about 1.2 and 1.5 
respectively. The 2 M & constraint is more restrictive. In this case the Hs, 
Qs, Qm, and Qh EOS are ruled out, while for Hm and Hh C14 must be less 
than about 0.55 and 1.16 respectively. 
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As the maximum observed neutron star mass is pushed upwards, and 
more is learned about the nuclear EOS, the observational upper bound on 
C14 will come down. If we assume the existence of a non-rotating neutron 
star of 2M , then even for the hardest equation of state we have considered 
we obtain the bound cu < 1.16. 

4.2 Surface redshift constraints 

There is not yet a definitive detection of an atomic spectral line from the 
surface of a neutron star. The strongest current case comes from stacked 
observations of thermonuclear X-ray bursts from EXO 0748-676, from which 
a surface redshift of 0.35 was inferred jlS] based on identification of some 
absorption-like features as being produced by highly ionized iron. The mass 
of this star is not certain, but Ozel [49] used simplifying assumptions about 
the constancy of the peak flux of the bursts and radiative transfer to suggest 
that the mass of this object is probably not less than 1.8 M . 

Measurements such as these, once confirmed, can provide a joint con- 
straint on C14 and the equation of state via the dependence of surface red- 
shift on mass in ae-theory. (We have checked that the mass inferred using the 
method of Ref. [49j differs from the GR value by less than 2% when cu = 1, 
so the leading order effect of cu is only in the relation between radius and 
mass, equivalently redshift and mass.) Fig. 0] shows a plot of z versus c u for 
1.8 Mq stars using the Hm, Hh, and Qh EOS. These are the three hardest 
equation of state models and have equilibrium configurations at this mass. 
(As shown in Fig. [31 the other three softest equation of state models do not 
have equilibrium configurations at this mass.) The surface redshifts increase 
by roughly 10% as C14 ranges from to 1. 

If in the future surface redshifts together with masses can reliably be de- 
termined, then tight constraints on the equation of state in GR may be ob- 
tained by combining measurements for a collection of stars. It is also possible 
that single measurements may provide stringent constraints. For example, as 
revealed in Fig. HJ the proposed surface redshift 0.35 of EXO 0748-676 [4S] 
is compatible with 1.8 M Q only for the hardest eos (Hh) among those we 
considered. In general, the parameter cu could not be constrained without 
separate knowledge of the equation of state. However, in the example just 
mentioned one could serendipitously tightly constrain both the equation of 
state and the value of C14, since a redshift of 0.35 is the lower limit of all the 
curves at this mass. 
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Figure 4: Surface redshift factor z versus Cu for 1.8 M Q neutron stars using 
the hardest equations of state. Hm (solid) is on top, Qh (dashed) is in the 
middle and Hh (solid) is on the bottom. Note that the GR value of z = 0.35 
for the hardest EOS, Hh, is consistent with the proposed surface redshift 
of 0.35 for EXO 0748-676 [48]. The Hm and Qh lines begin to curve up 
near Cx4 =1.1-1.2 because the maximum mass for these equations of state is 
approaching 1.8 M . 
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5 Discussion 



We have found that the structure of non-rotating neutron stars in Einstein- 
aether theory is fairly close to that in general relativity, but there are quanti- 
tative differences. Depending on the equation of state, the maximum masses 
range from about 6-15% smaller than in GR when the ae-theory parameter 
C14 is equal to 1. The corresponding surface redshifts are roughly 10% larger 
than in GR. Measurements of high gravitational masses or precise surface 
redshifts have the potential to yield strong joint constraints on the equation 
of state and on deviations from GR. Therefore, as laboratory experiments 
and other observations narrow down the equation of state of cold matter at 
several times nuclear density, neutron star observations may be a valuable 
resource for exploring deviations from general relativity in strong gravity. 

We now make some comments about related further work that would be 
interesting to pursue. First, the neutron star solutions we considered are 
at rest with respect to the asymptotic aether. Although corrections due to 
motion with respect to the aether are not significant for the present paper, 
they are important for the high precision predictions of radiation damping 
in compact binaries. In particular, the missing ingredient in the analysis 
of [3l] is the value of the "sensitivity" parameter measuring the velocity 
dependence of the mass. It should be possible to compute this parameter for 
different masses and different equations of state by determining the velocity 
perturbations of the solutions found here, or by finding the exact nonlinear 
solutions with finite velocity. 

We considered in this paper only neutron star phenomenology. What is 
the situation for black holes? It was found in [22] that although different 
from stellar exteriors, nonrotating black hole solutions in ae-theory are very 
similar to the Schwarzschild solution of GR. Hence it is not likely that these 
could lead to significant constraints. 

Nevertheless, it is quite conceivable that strong deviations from GR will 
be found for rapidly rotating black hole solutions. This is suggested by the 
presence of the ergoregion, in which the inertial frames are strongly dragged. 
The preferred frame aspects of ae-theory may be conspicuous here, due to 
larger gradients in the aether field. Also, unlike for the spherically symmetric 
non-rotating case, the spin-1 degrees of freedom of the aether could be acti- 
vated in the axially symmetric setting. To explore these issues would require 
finding numerical solutions describing rotating black holes in Einstein-aether 
theory. 
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